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Abstract 

We evaluate the antiproton-to-proton flux ratio which can be generated by 
neutralino-neutralino annihilation in the galactic halo, considering the most general com- 
positions for the relic neutralinos and modelling the neutralino local density according 
to its relic abundance. We find that in the case of mixed neutralino compositions this 
p/p ratio is higher than for pure higgsinos or gauginos. It is shown how the expected 
improvements in sensitivity for the new measurements of the p/p may provide very useful 
information, complementary to the one obtainable with other experimental means. 
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1. Introduction 



The neutralino (x) is generally considered the most favorite candidate for cold dark 
matter, and many possible ways for finding some evidence for it as a relic particle in the 
Universe have been envisaged. The most straightforward mean is obviously provided by the 
direct method, where the incoming of a dark matter (DM) neutralino would be signalled 
by the energy released, in an appropriate detector, by this particle when it scatters off a 
nucleus of the experimental apparatus [1-2]. The existence of relic neutralinos could also 
be inferred indirectly, by detecting some specific signals due to neutralino-neutralino anni- 
hilations. The best chances of obtaining indirect information of this sort would be provided 
by neutrino telescopes aiming at the measurement of fluxes of high-energy neutrinos from 
the Sun or from the centre of the Earth [3-4]. In fact these macroscopic bodies could effi- 
ciently capture and accumulate relic neutralinos, and finally emit the neutrinos produced 
by the neutralino pair annihilation in their interior. To be sensitive to this measurement 
the neutrino telescopes have to be of very large area (> 10 5 m 2 )[4]. Other indirect meth- 
ods consist in the detection of signals due to x~X annihilation occurring in our galactic 
halo [5] . Theoretical predictions for the magnitude of these signals have been performed in 
the case of the gamma-line [6-8] , of a diffuse gamma-ray flux [7-9] , of a gamma-ray flux 
from the center of the Galaxy [10,8] and of p, e + components in the cosmic rays [11-14]. 

The issue at stake, i.e. the detection of a DM particle in our Universe, is so important 
that every method of detection is worth to be pursued. This is even more true, since 
the odds of the various experimental means depend on the actual nature (composition 
and mass) of the neutralino. In the present paper we concentrate on the evaluation of the 
amount of antiprotons in cosmic rays due to the x~X annihilation in our halo. This problem 
has been previously addressed by other authors only for neutralinos of pure compositions: 
either higgsinos or gauginos, and for assigned values of the neutralino relic abundance [11- 
14] . At variance with these previous investigations we consider here neutralinos of the most 
general compositions and we model the neutralino local density according to its calculated 
relic abundance. Our results show how sensitive may be the antiproton-to-proton flux 
ratio, p/p, to the nature of the neutralino. We find that the experimental capabilities are 
not far from the possibility of a fruitful search for neutralino dark matter by measuring the 
p/p ratio. The chances of obtaining positive results in these experiments depend on the 
actual nature of the neutralino, these chances being higher for a x of mixed composition 
than for the exactly pure configurations. 

2 The p flux 

The differential rate (per unit volume) for production of p's from x~X annihilation is 
given by 



where Ep is the p energy, a is the x~X annihilation cross section, v the neutralino velocity 
in the halo, and p x is an averaged neutralino density to be discussed later; f(Ep) denotes 
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where F denotes the x~X annihilation final states, B^J is the branching ratio into quarks 

X.J 



or gluons in the channel F; dNp/dEp is the differential distribution of the antiprotons 
generated by hadronization of quarks (with the exception of the top quark) and of gluons. 

From Eq.(l) the standard formula for the interstellar p flux, due to x~X annihilation, 
follows 



where Vp and Tp are the p velocity and confinement time, respectively. 

Of the various quantities that enter in Eq.(3) the first two, < av > and f(Ep), 
depend on the neutralino properties, whereas the last two, Vp and Tp, are only related to 
the propagation of p's inside the galactic halo. The quantity p x depends on the intrinsic 
properties of x, but also on the features of the p propagation, since the average has to be 
taken over the region where the p's diffuse. In what follows we illustrate how we evaluated 
all these quantities. 

2.1 The x~ X annihilation cross section 

In the previous works [11-14] < av > has been estimated by assigning an arbitrary 
value to the neutralino relic abundance O x /i 2 and by exploiting a standard expression 
relating the relic abundance to the x~X annihilation cross section integrated from the 
freeze-out temperature to the present temperature [15]. Instead, we calculated < av > in 
the framework of a particle model, the Minimal Supersymmetric Standard Model (MSSM). 
For the evaluation of a we have included all the relevant diagrams at the tree-level, i.e. 
Higgs-exchange diagrams and Z-exchange diagrams in the s-channel and diagrams with 
exchanges of squarks, neutralinos and charginos in the t-channel, with all the final states 
which may contribute to the production of p's. We have also included the one-loop dia- 
grams which generate a two-gluon final state [16], since, as shown in Ref.[14], this con- 
tribution is important for p production by x~X annihilation in case of neutralinos with 
gaugino-dominance and with a mass m x lighter than the top mass. 

Our calculations have been performed for neutralinos of the most general composition 
in terms of the photino, zino and higgsino fields 
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7 and Z are related to the U(l) and SU(2) neutral gauginos, B and W 3 , by the standard 
rotation in terms of the Weinberg angle. For other details about the theoretical framework 
adopted here for our calculations we refer to our previous paper [17]. We only remind that 
the neutralino sector is usually described in terms of three parameters: M2, \i and tan/?. 
In the present paper for tan/3 we consider two representative values tan/3= 2, 8 whereas 
the other two parameters are varied in the ranges: 10 GeV < M 2 < 6 TeV, 10 GeV < \fi\ < 
3 TeV. 

We introduce a composition parameter P = a\ 2 + CI2 2 which gives the gaugino frac- 
tional weight. Almost everywhere in the parameter space, at fixed tan/3, there is a one- 
to-one correspondence: (M 2 ,/u) <-> (P,m x ); thus in what follows we replace the two pa- 
rameters (M2, fS) with (P, m x ), since these variables have a more direct physical meaning 
for the neutralino. 

As for the unknown masses (i.e. masses for Higgs bosons, Susy particles and top 
quark), it has to be specified that most of the results of our analysis presented here refer 
to the following set of values: m t = 170 GeV, = 50 GeV, mj = 1.2 m x , when m x > 45 
GeV, rrij = 45 GeV otherwise, except for the mass of the top scalar partners (the only 
ones relevant to radiative corrections) which has been taken m = 1 TeV. Only in a few 
cases, other choices, to be defined later, have been adopted for some of these masses. 

2.2 The p differential spectrum 

Let us consider now the p differential distribution f(Ep). The branching ratios B^p 
have been calculated for all annihilation final states which may produce p's. These states 
fall into two categories: i) direct production of quarks and gluons, ii) generation of quarks 
through intermediate production of Higgs bosons, gauge bosons and t quark. In order to 
obtain the distributions dNp /dEp the hadronization of quarks (with the exception of the 
top quark) and gluons has been evaluated by using the Monte Carlo code Jetset 7.2 [18,19]. 

The ensuing differential distributions show the usual scaling properties in the variable 
x = Ep/(E cm /2), where E cm is the center of mass energy [20,13]. However, sizeable 
violations of this scaling behavior occur for small values of x (x < 0.1) [20]. Since this 
is precisely the energy range of interest for the process under discussion, instead of using 
these scaling properties, we have evaluated the differential distributions at four different 
values: E cm /2 = 15, 30, 60, 120 GeV. The differential distributions at different E cm , have 
been obtained by interpolation. A convenient analytic expression that fits our numerical 
spectra is provided by the following form 
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where w = Log(x c ), x c = Tp/(E cm /2) and Tp is the p kinetic energy. The parameters p^'s 
are given in Table 1 for the value (E cm /2) = 30 GeV. 

As an example, our results for the case of the differential distribution due to the 
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b quark fragmentation are shown in Fig.l in terms of the variable x c . Our spectra are 
numerically in agreement with the ones of Ref.[13] for the cases considered in that paper, 
i.e. for the fragmentations of the c and b quarks. 

2.3 Propagation of the p's in the halo 

Let us turn now to the evaluation of the quantities that depend on properties of the p 
propagation in the halo. In the diffusion model [21] for propagation of cosmic rays in our 
Galaxy the confinement time Tp may be evaluated from the formula 

_ hl l-\{h g /h h f + \{h g /h h f ^ hi 
P 3D l-\{h g /h h ) ~3D W 

where h g and hh are the half-thicknesses of the disk and of the confinement region, re- 
spectively, and D is the diffusion coefficient. The approximation involved in Eq.(6) is 
justified since h g = 100 — 150 pc and hh > 700 pc. By varying hh over the range 
1 Kpc < hh < 10 Kpc and by using for D the numerical values of Ref.[21] (the ratio 
D/hh is rather stable with values in the range D/hh — (1-7 — 3.3) • 10 6 cm s _1 ), one 
obtains for the p confinement time 

0.3 • 10 15 s < Tp < 6 • 10 15 s (7) 

in agreement with previous evaluations [12,13]. In the following, for definiteness, we adopt 
the value t p = 3 • 10 15 s. 

Let us discuss now the quantity p x appearing in Eq.(3). This represents the average 
value of the neutralino density over a sphere centered at the Earth position, with a radius 
of the order of hh- It is then straightforward to relate p x to the local (solar neighborhood) 
neutralino density p x ,ioc- Assuming for the DM halo the usual distribution p(r) = pioc(l + 
rQ/a 2 )/(l + r 2 /a 2 ), where r & ~ 8.5 Kpc and p\ oc ~ 0.3 GeV cm -3 one finds in general 

(< p 2 > /pL) V2 - 1-8 for a = 2.0 Kpc [22] (8) 

(< p 2 > /p 2 oc ) 1/2 ^ 0.97 for a = 7.8 Kpc [23]. (9) 

Here, to be conservative, we adopt an estimate close to the one in Eq.(9), p ~ p\ oc , which 
entails a similar relationship for the fraction of DM due to neutralinos, i.e. p x ~ P x ,ioc- 

2.4 Evaluation of p Xy \ oc 

In order to assign a value to p x ,\ oc , for any point of the parameter space (specified by a 
set of values for M 2 , p and tan/3), we have evaluated the neutralino relic abundance Q x h 2 
and we have then determined the value R = min{l, (Q x /i 2 )/(Q/i 2 ) min }, where (fi/i 2 ) min is 
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defined as the minimum value of (fi/i 2 ) compatible with observations. Typically, taking 
> 0.1 and h ~ 0.5, one has (0/i 2 ) min ~ 0.03; this is in fact the value employed here for 
(fi/i 2 ) min . Then p x ,i oc has been determined by the standard rescaling procedure 



P X ,loc = RP\oc • (10) 

The method adopted for the evaluation of Q x h 2 is the one illustrated in Ref.[17]. We 
only remind here that O x /i 2 is inversely proportional to the averaged annihilation cross 
section < av >, integrated from the freeze-out temperature to the present one. Thus it 
follows that when, in a given point of the parameter space, a is large, in that point Q x h 2 
is depleted. 

2.5 Modulation of the p flux 

The p flux, measured at Earth, may be obtained from the interstellar p flux given 
by Eq.(3) by taking into account the modulation due to the solar wind. This effect is 
calculated here by employing the procedure discussed in Ref.s [24,25] and also utilized 
in Ref.s [11-14]. Denoting by E and T the p total and kinetic energies and by k the p 
momentum measurable at Earth, the solar-modulated p flux is given by 

**W = W^^viEp) . (11) 

E and Ep are related by 



k + E 

E p = k c In +E C + AE 

Ev = E + AE 



for k < k c 
for k > k r 



(12) 



The numerical values for k c and AE depend on the phase of the 11 year solar cycle. In 
what follows only the antiproton-to-proton flux ratio p/p, rather than the individual p 
and p fluxes, will be considered. Since these two fluxes are modulated in the same way, p/p 
is quite insensitive to the specific values adopted for k c and AE. To be definite we use in 
the following the values corresponding to a period of minimum solar activity: k c = 1.015 
GeV and AE = 495 MeV. 

3 Results and conclusions 

From Eq.s(3) and (10) we see that, at fixed m x , the flux $p depends on the neutralino 
properties mainly through the quantity < av > R 2 (some further, less effective dependence, 
is introduced in $p through f(Ep)). Thus we start the presentation of our results by 
discussing the behavior of < av > R 2 as we move in the neutralino parameter space. In 
Fig. 2 we display < av >, R 2 and their product < av > R 2 as functions of P, for m x = 30 
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GeV, tan/? = 2, 8 and for the set of mass parameters (m^ and mj) previously stated. In 
this figure we note a number of interesting properties: i) the x~X annihilation cross section 
is maximal for large higgsino-gaugino mixing (P ~ 0.5) (this feature is expected, since 
the small value used here for rrih makes the Higgs-exchange amplitude of the annihilation 
cross section very large and dominant over the other amplitudes; this effect is particularly 
important when x is largely mixed); ii) the neutralino relic abundance fi x /i 2 , and then 
also R 2 , are strongly depleted except for rather pure compositions, since the annihilation 
cross section is very efficient unless P ~ or P ~ 1 (for the reason mentioned at point 
i)); iii) from the plot of < av > R 2 one sees that the combination of the two previous 
properties makes the neutralinos with a mixed composition ( P ~ 10~ 3 — 10~ 2 in the 
present representative points) as the most favourite ones for detection by the method 
under discussion; iv) exactly pure compositions (P = 0,1) provide rather low values for 

< av > R 2 and are in general not indicative of what the real situation could be when some 
mixing is turned on. 

In Fig. 3 we show how the quantity < av > R 2 depends on the values assigned to 
the masses rrih and mj. The dashed lines correspond to the values: i) rrih = 50 GeV, 
raj = 150 GeV in Fig.3a, ii) m h = 120 GeV, mj = 45 GeV in Fig.3b, iii) m h = 120 
GeV, rrij = 150 GeV in Fig. 3c. The solid line always denotes the case where rrih = 50 
GeV, mj = 45 GeV; this line, already shown in Fig. 2, is displayed here only for sake 
of comparison. In Fig. 3a we see that a marked difference between the two curves shows 
up at a mixed composition with gaugino dominance. This is expected since in this case 
sfermion-exchange diagrams are important; then, for the gaugino-dominated neutralinos, 
an increase of m i entails a sizeable suppression in the x~X annihilation cross section with a 
subsequent enhancement of the relic abundance and then of the R 2 scaling factor. Fig. 3b 
shows a similar effect of enhancement; however, here the increase in R 2 , and then in 

< av > R 2 , occurs for the mixed compositions, since these are very sensitive to the value 
employed for rrih- Finally in Fig. 3c we notice that, when rrih and rrij are both large, the 
curve has a marked peak for P ~ 0.5; as a matter of fact in this case rescaling is practically 
absent and < av > R 2 ~< av >. 

In Fig. 4 the ratio of antiproton-to-proton integrated fluxes, (p/p)int, is displayed 
versus P for m x = 30 GeV and tan (3 = 2,8. The p spectrum has been eval- 
uated from Eq.(3) and the proton spectrum is taken from Ref. [26]: <& p (E p ) = 
1.93 v p (E p /GeV)~ 2 ' 7 cm _2 s _1 sr _1 GeV _1 (v p and E p are the proton velocity and en- 
ergy). The modulation for both spectra is obtained with the procedure discussed above in 
Sect. 2. 5. Both p and p spectra have been integrated over the ranges: (a) 100 MeV< T < 
640 MeV; (b) 640 MeV< T < 1580 MeV; and (c) 120 MeV< T <860 MeV in order to con- 
form to the energy windows used in the experiments of Ref.s [27,28]. The horizontal lines 
in Fig. 4 denote the corresponding upper bounds: (p/fOint = 2.8 • 10 -5 (85 % C.L.) for the 
interval a) [27]; (p/p) int = 6.1 -lO" 5 (85 % C.L.) for the interval b) [27]; (p/p) int = 1.8- 10" 5 
(85 % C.L.) for the interval c) [28]. It is remarkable that the experimental sensitivities are 
already at the level of the maximum of the theoretical curves (at the representative points 
under discussion). This means that by improving these sensitivities one is able to start 
investigating the relic neutralinos by the measurement of p/p. 
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Let us turn to the plots of the p/p ratio versus the p kinetic energy T. Fig. 5 shows how 
sensitively, at fixed m x , this ratio depends on the neutralino composition. In this figure 
also the background due to the secondary antiprotons produced in cosmic-ray collisions in 
the interstellar medium is displayed [29]. 

The dependence on the neutralino mass of the plots of p/p versus T is illustrated in 
Fig. 6 for three values of P separately [30]. In Fig. 6c we notice that the ordering of the four 
curves for our signal is the usual one, expected on the basis of the dependence on m~ 2 in 
Eq.(3). However, in Fig.s 6a,6b this mass hierarchy turns out to be modified; this effect 
occurs when 2m x is close to one of the masses of the particles exchanged in the s-channel 
of the annihilation process. Then, in the present case, this happens for m x = 30 GeV 
which is close to the value used here for one of the neutral Higgs bosons. 

It has to be pointed out that many uncertainties are involved in the theoretical pre- 
dictions for the p/p signal, as it emerges from the derivation of the p flux presented above. 
Suffice it to mention the number of free parameters in MSSM and the wide range in the 
estimate of the confinement time Tp. 

Nevertheless, from the previous analysis it turns out that a substantial improvement 
in the measurement of the p/p ratio [31] can make this kind of investigation competitive 
in the future with other experimental strategies which aim at finding a sign of neutralino 
DM. This measurement shares with the measurements of other signals due to the x~ 
X annihilation in the halo [6-10] the peculiar possibility of providing information about 
^-compositions with small mixings, which have little chances to be detected by other 
experimental means [32]. This point is illustrated in Fig. 7 where we give the (p/p)i n t ratio 
versus m x (7a) and versus P (7b) in the form of a scatter plot, obtained by varying the 
two parameters M2 and \x in the ranges 10 GeV < M2 < 6 TeV, 10 GeV < < 3 
TeV, at fixed tan/3 = 8. These two plots show how the exploration may proceed inside 
the DM neutralino physical region, especially at small P and small m x , by increasing the 
experimental sensitivity in the measurement of p/p. 

The most serious problem for the detection of the neutralino signal under discussion 
is due to the uncertainty in the evaluation of the flux of secondary p's produced by cosmic- 
ray collisions; this flux plays here the role of a background. Unfortunately, the calculation 
of this spectrum is affected by large uncertainties at small p kinetic energies, making the 
signal-to-background discrimination difficult. Then a more accurate evaluation of this 
background is very much needed. However, it has to be emphasized that new stringent 
upper limits in the range of small energies would anyway provide useful constraints about 
the neutralino physical region. 
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Figure Captions 

Figure 1 p differential distibution dNp / dx c from the fragmentation of a b quark, as 
a function of x c = Tp/(E cm /2), Tp being the antiproton kinetic energy. The curves refer 
to four different values of E cra . From top to bottom: E cul /2 = 120, 60, 30 and 15 GeV. 

Figure 2 < av >, R 2 and < av > R 2 as functions of the neutralino composition at 
a fixed value of the neutralino mass m x = 30 GeV. The other parameters are: tan (3 = 8 
(solid curve), tan/3 = 2 (dashed curve), = 50 GeV, mj = 45 GeV. In these plots, in 
order to represent the dependence of various quantities accurately for very small mixings 
(i.e. for P ~ 0, P ~ 1), on the left part of the horizontal axis the variable P is given, 
whereas on the right part of the same axis the complementary variable 1 — P is reported, 
both in a log-scale. 

Figure 3 < av > R 2 as a function of the neutralino composition for tan/3 = 8 and 
for a fixed value for the neutralino mass m x = 30 GeV. The solid curve refers to rrih = 50 
GeV, mj = 45 GeV. The dashed curves denote different choices for the mass parameters: 
(a) m h = 50 GeV and m f ~ = 150 GeV; (b) m h = 120 GeV and m f ~ = 45 GeV; (c) m h = 120 
GeV and m~ f = 150 GeV. 

Figure 4 Ratio of antiproton-to-proton integrated fluxes, (p/p)i n t, as a function of 
the neutralino composition for m x = 30 GeV. The solid curve refers to tan/3 = 8, the 
dashed curve to tan/3 = 2. The mass parameters are: mu = 50 GeV and rnj = 45 GeV. 
The integration ranges are: (a) 100 MeV< T < 640 MeV; (b) 640 MeV< T < 1580 MeV; 
(c) 120 MeV< T <860 MeV. The horizontal lines are the corrisponding present upper 
bounds: (a) (p/p) int = 2.8 • 10" 5 (85 % C.L.) [27]; (b) (p/p) int = 6.1 • 10" 5 (85 % C.L.) 
[27]; (c) (p/p) int = 1.8 • lO" 5 (85 % C.L.) [28]. 

Figure 5 Antiproton-to-proton flux ratio, p/p, as a function of the kinetic energy T, 
for m x = 30 GeV and for different neutralino compositions: P = 10~ 3 (long-short dashed 
curve); P = 0.1 (solid curve); P = 0.5 (long dashed curve); P = 0.9 (dash-dotted curve); 
P = 0.999 (short dashed curve). The dotted curve is the background flux of secondary 
antiprotons, calculated in the leaky-box model [29]. Here tan/3 = 8, mh = 50 GeV and 
mj = 45GeV. 

Figure 6 Antiproton-to-proton flux ratio, p/p, as a function of the kinetic energy 
T. Different curves refer to the following neutralino masses: m x = 15 GeV (solid curve); 
m x = 30 GeV (dashed curve); m x = 60 GeV (dash-dotted curve); m x = 120 GeV (long- 
short dashed curve). In (a) the neutralino composition is fixed at the value P = 10 -3 ; 
in (b) P = 0.5; in (c) P = 0.999. The dotted curve is the background flux of secondary 
antiprotons, calculated in the leaky-box model [29]. Here tan/3 = 8, rrih = 50 GeV and 
mj = 1.2 m x , when m x > 45 GeV, mj = 45 GeV otherwise. 

Figure 7 Ratio of antiproton-to-proton integrated fluxes, (p/p)\ n u as a function of 
the neutralino mass m x (a) and of the neutralino composition (b). The figure is presented 
in the form of a scatter plot obtained by varying the model parameters in the ranges: 10 
GeV < M 2 < 6 TeV, 10 GeV < \p\ < 3 TeV. The integration interval is 120 MeV< T <860 
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MeV. The horizontal line is the corresponding upper bound (p/p)i n t = 1-8 • 10~ 5 (85 % 
C.L.) [28]. Here tan/3 = 8, = 50 GeV and mj = 1.2 m x , when m x > 45 GeV, mj = 45 
GeV otherwise. 
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Table Caption 



Tablel Values for the p^'s parameters of the distribution given in Eq.(5) for the value 
-Ecm/2 = 30 GeV. Each row corresponds to a quark flavour and to the gluon. 





Pi 


P2 


P3 


P4 


u, d 


1.7443 


2.9722 


-0.62737 


-1.1193 


s 


7.5915 


4.0079 


1.5257 


-1.4500 


c 


36.268 


4.8055 


3.0725 


-1.5923 


b 


6.9350 


5.4031 


1.7076 


-1.3286 


gluon 


132.86 


4.9732 


3.5982 


-1.6360 
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